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Half-Power Bandwidth method for the 
evaluation of synchronous and 


Nomenclature 

Ava Rotor synchronous response amplitude and 
phase, respectively 

Face Amplitude A at resonance 

BB Rotor amplitude and phase in response to 
nonsynchronous excitation 

Co Rotor 2x response amplitude and phase, 
respectively 

D,K,M Rotor first lateral mode modal damping, stiff- 
ness and mass, respectively 

e = 2.7183 

E,8 Cracked rotor synchronous response vector 
components 

F.y Nonsynchronous perturbation force ampli- 
tude and phase, respectively 

j =¥-1 

m ré Rotor modal unbalance mass, radius, and 
angular orientation 

P,é Constant radial force amplitude and angular 
orientation 

2 OP Nonsynchronous amplification factor 

Qum Synchronous amplification factor 

i RER Synchronous amplification factor for a 
cracked shaft 

Qa Amplification factor for 2x response 
component 

A Difference of rotor stiffness in two orthogonal 
directions perpendicular to rotor axis, a mea- 
sure of a rotor crack 

Č Modal damping factor 

A Fluid average circumferential velocity ratio 

Ww External perturbation frequency 

Wes Resonance frequency 

Wo, 155 Frequencies at the rotor response phase lag of 
45° and 135°, respectively 

Q Rotative speed 

Raa VKM First balance resonance speed 

Ql oa Rotative speed at 1X resonance of an 
anisotropic rotor 

E Rotative speed at resonance of 2X component 
axosbt a5 Rotative speeds at the 1X response phase lag 
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ghis article discusses the use of the Half-Power Band- 
width (HPB) method for the evaluation of the syn- 
chronous and nonsynchronous quadrature dynamic 
stiffnesses and amplification factors for rotor/bearing/seal sys- 
tems. Traditionally, the HPB method is referred to as the 
modal damping identification method. This article shows that 
this method not only provides damping, but also the entire 
quadrature dynamic stiffness, which, in the case of rotor lateral 
modes, also contains tangential, cross-coupled stiffness terms. 
The HPB method, when used on rotor nonsynchronous per- 
turbation test data, vields an evaluation of the destabilizing, 
flow-induced forces resulting from the fluid environment. 
When the HPB method is used for rotor 2X response, it may 
provide information of an impending shaft crack. 


Synchronous response of an isotropic rotor in the 
range of its first lateral mode 

As has been shown, [References 2-5], the synchronous 
response vector, due to unbalance of a simple isotropic rotor, 
for either vertical or horizontal motion in the range of the 
rotor first lateral mode, is as follows: 
7 mri*e!* 

K-MQ* + jQD(1-A) 

where Ae" is the 1X response vector (A = amplitude, a = 
phase), K, M, Dare modal stiffness, mass, and damping; m, r, & 


Ae! j= V (1) 
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are modal unbalance mass (“heavy spot’), radius, and angular 
orientation, respectively, & is the rotative speed and A is the 
fluid average circumferential velocity ratio, the cross-coupled 
tangential effect of the fluid surrounding the rotor [2]. The 
product D(1-A), is referred to as “modal effective damping.” 


In Equation (1) the expression in the denominator repre- 
sents the system complex dynamic stiffness with its direct 
term, K-Mf", and quadrature term, 2D(1-A). 

In the “real” (direct) or “imaginary” (quadrature) plane, 
the response vector (1) represents a point for any constant 
rotative speed, @. For the speeds varying from zero to some 
value higher than the resonant one, the response vector draws 
a figure similar to a circle (Figure 1), known as a modal circle, 
“Nyquist circle,” or more commonly, a “polar plot.” 


There are several characteristic points on the polar plot. At 
Q=0, there is A = 0 and a= 4; thus, the polar plot begins at 
the coordinate origin with the response vector’s phase, 6, the 
same as the heavy spot phase. For low values of damping, at 
Q= Q, = VKM, the response vector reaches its highest, reso- 
nant value, and the phase lags the original value 6(at Q = 0), by 
90 degrees. At high rotative speed, the inertia term, MQ*, in 
Equation(]) becomes dominant, and the response vector 
amplitude tends to a constant value, while the phase lags 180 
degrees from its value at Q = 0. 


There are two more characteristic points on the polar plot, 
namely at a= 6-45° and a= 6-135". These points are situ- 
ated on the polar plot circle diameter which is perpendicular 
to the response vector at resonance (Figure 1), 

The values of the rotative speed at these points, as well as at 
resonance, leads to the basic equation of the half-power band- 
width method: 


Q 135 ilas nR DiI A) 
Qes VRM 
where €=D/(2@VKM) is known as a “damping factor.” The 


inverse of expression (2) represents the basic equation for the 
synchronous amplification factor Q 


0 eee 1 
ilias - ils iss 2€(1-A) 





= 2Z(1-A) (2) 


“yn! 


(3) 


Thus, in order to evaluate the rotor first mode effective 
damping and corresponding synchronous amplification fac- 
tor, three speed values from the polar plot are required: the 
resonant speed at the highest value of the response amplitude 
and two speeds at the polar plot circle diameter perpendicular 
to the response vector at resonance (Figure 1). The calculation 
is performed using Equations (2) or (3). 
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Trending the synchronous amplification factor may provide 
valuable information about possible changes in the damping 
factor or in fluid-induced forces, represented by the ratio A, 
which, if increased to one, may lead to rotor instability. Note 
that the synchronous amplification factor does not depend on 
the exciting unbalance force, provided that the response 
amplitude stays within a linear response range. 

The amplification factor can also be estimated from the 1X 
response Bode plot as a ratio of the response amplitude at 
resonance to the amplitude at high rotative speed (lower, 
however, than the second mode frequency range). This eval- 
uation is, however, less accurate. 


Short history of “Half-Power Bandwidth" 

The name “Half-Power Bandwidth” originated from the 
response vector voltage (V) of an electrical circuit, which is 
usually plotted in logarithmic scale. The response amplitudes 
at a= $6-45°%and a=6-135° have approximately the same 
value, which is a 12 fraction of the resonant amplitude: 

a . 
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Figure 1 
Polar plot of the synchronous (1X) response of a rotor, 
measured by a lateral probe. 


June 1994 


Since power is proportional to the square of the voltage, the 
As- and A 5- points also represent the half power points of 
voltage. 

The voltage or displacement amplitude ratio 1/72 corre- 
sponds to approximately a 70% or 3dB amplitude reduction. 
Thus, the points defined by that ratio are sometimes called the 
70% or 3dB amplitude points, as well as the half-power points. 
These points define the 3db amplitude bandwidth used to 
calculate Q, 


asym 

With current data acquisition technology, the points can be 
determined most accurately from the numerical listing of the 
data. The second choice for accurate data would be the polar 
plot, which typically offers better resolution of the phase and 
frequency data than the Bode plot. 


Nonsynchronous response of an isotropic rotor with 
fluid interaction in the range of its first lateral mode 

Nonsynchronous excitation of a rotor occurs when an 
external-perturbation, rotating force with frequency w (w#Q) 
is applied, while the rotor rotates at a constant rotative 
speed, Q. 


The response vector, Bef, within the first lateral mode of a 
very simple isotropic rotor excited by the nonsynchronous 
rotating force, is as follows [3]; 


| Fe” 
EE aa E (4) 
K-Mw? + jD(w-AQ) 


where F and y are exciting force amplitude and angular orien- 
tation, respectively. Now the cross-coupled, tangential effect 
of the fluid surrounding the rotor appears as a separate com- 
ponent, DAQ, of the dynamic stiffness. 

A resonance (“‘direct” or mechanical”) still exists at 
Wie, =VKIM. The modification occurs in the quadrature 
dynamic stiffness term. This term in Equation (4) can no 
longer be always positive, as it was in Equation (1). Due to the 
fluid interaction transferring rotative energy into the lateral 
vibration mode through the cross coupling stiffness term, 
DAQ, the entire quadrature dynamic stiffness may become 
zero. This leads to the quadrature resonance conditions at 
Woon 7 AQ in a simple rotor system [3,4]. 

Ifthe Half-Power Bandwidth analysis is applied to the rotor 
excited with nonsynchronous sweep frequency w, one can 
easily deduce that the nonsynchronous amplification factor 
for the direct resonance W ea = YKM is as follows: 


(ey = 


(5) 
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The half-power points correspond to the phases y—45° and 
y-135°. The nonsynchronous amplification factor, Quen: is 
magnified in comparison to Q, (Equation (3)), by the factor 
which depends on the closeness of the quadrature resonant 
frequency, Wq FAQ, to the direct resonant frequency, 
w,.,=VKIM . If both resonant frequencies coincide, that is, if 
Ag=VK0+2/M, then, obviously, an instability occurs, and 
the response amplitude B , as well as Q ayn: infinitely increase. 

From Equation (5), it is clearly seen that the Half-Power 
Bandwidth method applied to rotor responses provides more 
than a simple damping measure. It provides the measure of 
the entire quadrature dynamic stiffness, containing damping 
and the destabilizing, cross-coupled stiffness generated by the 
rotor rotational mode. For more complex systems, more 
terms may appear in the equation; however, the fundamental 
rules remain essentially the same. 


1X amplification factor of a cracked shaft 


The synchronous 1X response amplitude, E, and phase 9, of 
a transversely-cracked rotor excited by an unbalance force, can 
be presented as follows [5]; 





[(K-Ma*)cosé + Dosind}* + [(K-Mo*)siné-Dacoss]? 


(K-MQ*)(K-A-MQ*) + D709" 


E=mrn* 
(6) 


(K-Mo*)siné-Dacosé 
(K-4-M*)cosé + DQsiné 


where A is the difference between the rotor stiffnesses in two 
orthogonal directions perpendicular to the rotor axis. The 
parameter A can represent a measure of a shaft transversal 
crack. For clarity of presentation, the factor, A, has been 
omitted in Equations (6) and (7). As can be seen from Equa- 
tions (6) and (7), the cracked-shaft synchronous response 
depends on the location of the crack in relation to the heavy 
spot. The worst case results if they are separated by 45 
degrees, that is if 6 = 45°. This is the case that will be discussed 
in this article, 
The 1X resonance occurs when: 


| K-A/2 | 
= (8 
3 AA M ( ) 


which, for an impending crack, is slightly lower than the 
previously discussed resonant frequency, VKM. 

Using the Half-Power Bandwidth method, the synchronous 
amplification factor can be evaluated. Since at Q=0 the 
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assumed phase is @= 6 = 45°, the half-power frequency points 
to be calculated occur at @=0° and © =—90°, which are —45° 
and -135° from the origin. 


Using the Half-Power Bandwidth method, the synchronous 
amplification factor of a cracked rotor is as follows: 
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With an increasing value of A for an impending crack in the 
rotor, the amplification factor increases significantly. Its 
growth rate depends on the amount of damping in the system. 
The “crack ratio,” 4/2K, appears in the amplification factor as 
another “tangential component” subtracted from the damp- 
ing factor (Z), thus decreasing the denominator of the ampli- 
fication factor (9). It is necessary to point out that the growth 
of the amplification factor depends considerably on the crack 
versus unbalance location. Equation (9) defines the worst 
condition only. In any case, trending the synchronous ampli- 
fication factor may give information on an impending rotor 
crack. Note that, as before, the amplification factor does not 
depend on the unbalance force amplitude. 


2X amplification factor of a cracked shaft 

The 2X lateral response of the rotor may occur as the 
second harmonic of high 1X vibrations, due to nonlinearities 
in the system, The 2X response may, however, be entirely 
independent from 1X vibrations. It may be generated by a 
constant radial force applied to an anisotropic rotor. The rotor 
anisotropy may result from a rotor transversal crack. For the 
rotative speed, 9, the rotor 2X response vector, Ce, has, in 
this case, the following amplitude, C, and phase, g [5]: 


PA 


m | 
K-A 


C= — | 
(10) 
-2MQ* 





= 
2(2K-A) +R? 
DQ 


KIK-A' (11) 
= -2Mn* 
2K-A 


o =-é-arctan 





where A is a measure of a shaft crack, and P and € are constant 
radial force amplitude and angular orientation, respectively. 
The force, P, may result from rotor misalignment, from fluid 
flow side-load (as in single volute pumps) or from gravity (as 
on horizontal rotors). For clarity of presentation, the ratio A 
has been omitted in Equations (10) and (11). 

The 2X polar plot starts at @ = 0, with the phase o, =—e. At 
resonance, Le. when Q@=Q._.,=2VKM , the 2X response ampli- 
tude is the highest: C... = PMA/DK(K-A), with the phase 
O = -E-90°, At high rotative speeds, the amplitude 
approaches zero, with the phase approaching o,,=—eé-180°. 
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The 2X polar plot behaves, therefore, very similarly to the 1X 
polar plot circle. 


The application of the Half-Power Bandwidth method 
results in: 


Qiage — Rage D _ 


fl eed K(K-A) 7 
{aw 


Therefore, the amplification factor, Q,,, for the 2X response 


2(2K-A) 





<a 


af Ka 


Fa ZEKA eyn a3) 





This 2X amplification factor includes the effect of the rotor 
anisotropy, such as is generated by a shaft crack. For the 
isotropic rotor (when à = 0), there is Q,, = 22, which means 
Qa = nn (with A=0). If, for example, A=0.3K, then 
Qas = 0.910,,,,. If A = 0.5K, then Q,, =0.82Q,,,. While the 2X 
response amplitude (10) increases almost proportionally to 
the crack parameter A, surprisingly, the 2X amplification fac- 
tor does not increase, but decreases with A. It is also substan- 
tially lower than Qana. Trending the 2X amplification factor 
may, however, provide valuable information on an impending 
shaft crack. 


Conclusion 


The Half-Power Bandwidth method evaluates quadrature 
stiffness components. These may consist of destabilizing, 
fluid-induced, as well as shaft-crack-related, terms. Measure- 
ment and trending of nonsynchronous amplification factors 
on machines may provide information on changes in either 
fluid-induced factors or in shaft anisotropy. The latter, in turn, 
may be affected by shaft cracking. $ 
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